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Nomenclature 

A,  A. .  =  stretching  stiffnesses  for  transversely  isotropic  and  cross- 

■  J 

ply  orthotropic  plates 

a,b  =  plate  dimensions  in  x  and  y  directions 

B. B--  =  bending-stretching  coupling  stiffnesses  for  transversely 

•  J 

isotropic  and  cross-ply  orthotropic  plates 
Crs  =  coefficients  defined  in  Eqs.  (15) 

D.D^  =  bending  stiffnesses  for  transversely  isotropic  and  cross-ply 
orthotropic  plates 
dx  =  a(  )/3x 

Ec>Et  =  compressive  and  tensile  Young's  moduli  for  transversely 

isotropic  bimodulus  material 

E*,E  *  fiber  and  matrix  Young's  moduli 

f  m 

E1.E2.E3  =  Young's  moduli  in  directions  x,y,z 

F?  =  generalized  force  components  defined  in  Eqs.  (20) 

Gi3tG23  =  longitudinal-thickness  and  transverse-thickness  shear  moduli 
of  orthotropic  material 

G  ,G  .  =  thickness-shear  moduli  for  transversely  isotropic  bimodulus  material 

tv  £  L 

h  =  total  thickness  of  plate 

K  =  shear-correction  coefficient  for  transversely  isotropic  plate 


K4.K5 


=  shear-correction  coefficients  for  cross-ply  orthotropic  plate 


Mi*Ni 


matrix  coefficients  defined  in  Eqs.  (20) 

linear  differential  operators  defined  in  Eqs.  (8)1 

=  stress  couples  and  inplane  stress  resultants  /  J*lanoottHce<i 

J^Uncati0lL 

interpolation  function  at  node  i  |~L  - - — 


number  of  nodes  per  element 


£tjtrihutinH(/ 

JLyrit  1 nv i * « » _  _ 
1,*  /Avall*»d/or 


special 


mnwUnw  .lain 


Qc,Qt 


Q  .Q 

wx,yy 


^  i  j  k  a 


u.v.w 


=  compressive  and  tensile  plane-stress-reduced  stiffnesses  for 
isotropic  bimodulus  material 
=  thickness-shear  stress  resultants 

■  0/2)(Qc+Qt>.  Qc  -  Qt 

=  plane-stress-reduced  stiffnesses  for  orthotropic  bimodulus  material 
=  normal  pressure  and  its  peak  value 

=  thickness-shear  stiffnesses  for  transversely  isotropic  and  cross- 
ply  orthotropic  plates 

=  matrix  coefficients  defined  in  Eqs.  (20) 

=  midplane  displacement  coefficients  (amplitudes  of  u°,v°,w°) 


u.v.w  =  displacements  in  x,y,z  directions 


u°.v°.w°  =  midplane  displacements  in  x,y,z  directions 

Vf,v£  =  fiber  volume  fraction  and  effective  fiber  volume  fraction 

X.Y  =  bending-slope  coefficients  (amplitudes  of  ,\|>  ) 

x  y 

x,y,z  =  plate  coordinates  in  longitudinal,  transverse,  and  downward 


Z’Zx’Zy 


thickness  directions 

■  Vh>  znx/h-  Vh 

=  neutral-surface  position  for  isotropic  square  plate 


o 

ej’eJ 


z  ,z  =  neutral -surface  positions  associated  with  e  =0  and  e  =0 

MX  •■J'  X  Jr 

a, 6  =  7r/a,  u/b 

£.,£?  =  strain  component  at  arbitrary  location  and  at  midplane 

J  J 

k  •  =  curvature  component 

w 

vr.v  =  fiber  and  matrix  Poisson's  ratios 

f  m 

vi2»v23  =  major  (longitudinal-transverse)  and  transverse-thickness  Poisson's 

ratios 


stress  component 
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*  ,<t>^  =  typical  variable  in  general  and  its  value  at  node  i 

C  0 

<li  ,ip  -  bending  slopes  in  xz  and  yz  planes 
x  y 

(  ),x  ■  3(  )/3X 

Subscripts: 

i,j  =  1,2,6  contracted  indices 

k  =  l(t  or  tension),  2(c  or  compression) 

i  =  layer  number 

Introduction 

The  increasing  use  of  composite  materials  in  structures  has  led  to  the 
requirement  of  more  realistic  mathematical  modeling  of  the  material  behavior 
and  incorporation  of  this  more  realistic  model  into  structural  analyses.  It 
has  been  found  that  certain  fiber-reinforced  materials,  especially  those  with 
very  soft  matrices  (for  example,  cord-rubber  composites),  exhibit  quite  dif¬ 
ferent  elastic  behavior  depending  upon  whether  the  fiber-direction  strain  is 
tensile  or  compressive1"3.  As  a  first  approximation,  the  stress-strain  be¬ 
havior  of  such  materials  is  often  represented  as  being  bilinear,  with  differ¬ 
ent  slopes  (elastic  properties)  depending  upon  the  sign  of  the  fiber-direction 
strain.  Such  a  material  is  called  a  bimodulus  composite  material,  and  it  has 
been  shown  that  the  fiber-governed  symmetric-compliance  model  proposed  in 
Ref.  4  agrees  well  with  experimental  data  for  several  materials  with  drastically 
different  elastic  properties  in  tension  and  compression. 

To  the  best  of  the  present  investigators'  knowledge,  the  only  previous 
analyses  of  plates  laminated  of  bimodulus  composite  materials  are  all  limited 
to  thin  plates.  Jones  and  Morgan5  considered  cylindrical  bending  of  a  finite- 
width  cross-ply  strip;  Klncannon  et  al.6  considered  cross-ply  elliptic  plates. 
Rectangular  plates  were  treated  by  Bert  et  al.7  using  a  closed-form  solution 
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and  by  Reddy8  using  mixed  finite  elements. 

Apparently  the  only  previous  analyses  involving  thick  plates  of  bimodulus 
material  are  those  of  Shapiro9  using  a  stress-function  elasticity  approach  for 
isotropic  circular  plates  and  of  Kamiya10  using  an  energy  approach  for  cylindrical 
bending  of  finite-width  isotropic  strips. 

Of  course,  for  thick  plates  laminated  of  ordinary  (not  bimodulus)  materials, 
there  have  been  a  number  of  analyses,  such  as  those  of  Whitney11,  Whitney  and 
Pagano12,  and  Turvey13  ,  for  example. 

The  analyses  presented  here  are  believed  to  be  the  very  first  analyses  of 
thick  plates  that  are  finite  in  two  directions  and  laminated  of  bimodulus  com¬ 
posite  materials. 

Formulation 

The  basic  theory  of  laminated  anisotropic  thick  plates  used  by  Whitney 
and  Pagano12  is  an  extension  of  Reissner's  theory  for  isotropic  plates14.  It 
is  based  upon  the  following  assumed  displacement  field: 
u(x,y,z)  =  u°(x,y)  +  z  <Px(x,y) 

v(x,y,z)  =  v°(x,y)  +  z*  (x,y)  (1) 

w(x,y,z)  =  w°(x,y) 

Here  x,y  are  rectangular  coordinates  in  the  plane  of  the  plate,  z  is  the 

thickness-direction  coordinate  measured  downward  from  the  midplane  of  the 

plate;  u,v,w  are  the  displacements  in  the  respective  x,y,z  directions;  u°, 

v°,w°  are  the  corresponding  midplane  displacements;  and  *  and  p  are  the 

x  y 

slopes  in  the  xz  and  yz  planes  due  to  bending  only. 

Neglecting  body  forces,  body  moments,  and  surface  shearing  forces,  the 
equations  of  equilibrium  can  be  written  as 


1 


I 


1,  +  =  0  ; 

N 

+  N„ 

=  0 

i,x  6,y 

6  ,X 

2  ,y 

1  +  Q  +  q  =0 ; 

M 

+  M 

-  Q  = 

x,x  Ky,y  ^ 

6  ,X 

2,y 

y 

M  ,  +  M„  -  Q  =  0 

i,x  6,y  ^x 

Here  q  is  the  normal  pressure,  (  ),x  denotes  3(  )/3x, and 


(N.,M.)  = 


(Qy,Qx)  - 


.h/2 

-h/2 

h/2 

-h/2 


(l,z)o.  dz 


(o4,o5)  dz 


(1-1 .2,6) 


(2) 


(3) 


Here  h  is  the  plate  (laminate)  thickness, and  the  so-called  contracted  sub¬ 
script  notation  is  employed  to  denote  the  stress  components.  Thus,  aj  and 
02  are  inplane  normal  stresses  in  the  x  and  y  directions;  o6  is  the  in¬ 
plane  shear  stress  associated  with  the  x,y  axes;  and  o4  and  o5  are  the  thick¬ 
ness  shear  stresses  in  the  yz  and  xz  planes. 

Assuming  that  the  only  plane  of  symmetry  existing  is  in  the  plane  of  the 
plate,  the  plate  constitutive  relations  can  be  written  as 


N, 


1  J 


>  = 


A. .  B.. 

B.  .  D. . 

L  U  U 


(i,j=l,2,6) 


k4s44  k4k5s45 

K4K5S45  kISss 


r  \ 


e5 

V  J 


(4) 


The  A..,  B,.,  D..,  S. .  are  the  respective  inplane,  bending-inplane  coupling, 

I  J  I  J  I  J  I  J 

bending  or  twisting,  and  thickness-shear  stiffnesses  defined  as  follows: 

h/2 

(Aij’Bij’Dij*  =  L,  (1,z,z2)Qijkt  dz 


r 

■1. 


■h/2 

h/2 

h/2 


(5) 


Mjk*, 


dz  ( i ,j=4,5) 
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Here  Q.  .kA  denotes  the  plane-stress  reduced  stiffness,  where  ij  refers  to 
the  position  in  the  stress-strain-relation  array  (analogous  to  Eqs.  (4)),  k 
refers  to  the  sign  of  the  fiber-directicn  strain  (1  ^  +,  2  a.  -),  and  i  is 
the  layer  number. 

The  linear  equations  for  the  kinematics  of  deformation  are 


e°  =  u’x  » 

4  ■  '.J,  i 

4  ■  u.y  ♦  % 

=  ^x ,x  ; 

<2  "  *y,y  I 

'  *x,y  *  *y,x 

(6) 


£4  =  w,y  +  ;  e5  =  w,x  + 


Equations  (l)-(6)  plus  those  of  Appendix  A  constitute  the  appropriate 
theory,  in  differential-equation  form,  for  the  class  of  plates  considered 
here  (linear,  thick,  laminated,  anisotropic,  bimodulus). 

Closed-Form  Solution  for  Cross-Ply  Laminate 

Here  we  consider  the  particular  case  of  a  cross-ply  laminate,  i.e.,  one 

in  which  some  of  the  layers  are  oriented  along  the  x  axis  and  the  remainder 

along  the  y  axis.  Then  the  terms  with  subscripts  16,  26,  and  45  vanish  from 

the  symmetric  arrays  in  Eqs.  (4).  For  bimodulus-material  cross-ply  laminates, 

Eqs.  (5)  integrate  as  indicated  in  Appendix  A  and  depend  upon  the  neutral- 

surface  locations,  Z  and  Z  ,  as  well  as  the  Q..  . 

x  y  l  j  kx. 

If  it  is  tentatively  assumed  that  the  neutral -surface  locations  are  in¬ 
dependent  of  x  and  y,  Eqs.  (2),  (4),  and  (6)  can  be  combined  to  yield  the 
following  governing  equations  in  terms  of  the  midplane  displacements  (u°,v°,w°) 

and  bending  slopes  U  and  ) : 

y  x 

[Lrs]{u°,v°,w°,*y,*x}T  *  {0,0,q,0,0}T 


•J 


(7) 
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where  [L  ]  is  a  symmetric  linear  differential  operator  matrix  with  the  follow 
rs 

ing  elements: 

Ll 1  =  Alldx  +  A66dy  >  L12  =  (A12+A66)dxdy  i  L i 3  =  0  ; 

L 14  =  (Bi2+B66)dxdy  ;  L15  =  Bnd£  +  B66d^  ;  L22  =  A66dx  +  A22dy  * 

L23  =0  l  L24  =  B66d*  +  B22d^  ;  L25  =  L14  ;  L33  =  -  K^Sgsd^  - 
K4  S44dy  ;  L 34  =  -  K4S44dy  ;  L35  =  -  K3  S55dx  ;  ^ 

1*44  3  Bgedx  +  D22dy  ■  ^4^44  >  ^45  =  ( Dl  2+D6 6 ' dxdy  5 

1-55  =  ^lldx  +  ^6 6 dy  “  K5S55  ;  d^  =  3 (  )/3x  ;  dy  =  3(  )/3y 

For  a  plate  hinged  flexurally,  but  freetomove  in  a  direction  normal  to 
each  edge,  the  boundary  conditions  are 


M0,y)  =  Ma,y)  =  0 
v°(0,y)  =  v°(a,y)  =  0 
w°(0,y)  =  w°(a,y)  =  0 


u°(x,0)  =  u°(x,b)  =  0 


<f<y(0,y)  »  *  (a,y)  =  0 


N2(x,0)  =  N2(x,b)  =  0 

w°(x,0)  =  w°(x,b)  =  0  (9) 

M2(x,0)  =  M2(x,b)  =  0 

Mi(0,y)  =  Mi(a,y)  =  0  ; 

The  criteria  that  the  neutral -surface  locations  associated  with  the  x 
and  y  directions  remain  constant  are  as  follows: 

e?  +  znx<!  =  0  ;  c°2  +  znyK2  =  0  (10) 


*x(x,0)  =  4»x(x,b)  =  0 


or 


nx 


e°/<i  =  -  u°  /l|» 


Zny  =  -  *>2  -  ~  V> 


,x  Tx,x 

0  /* 

.y  vy,y 


(ll) 


i 

i 


! 


9 


The  normal-pressure  loading  is  taken  to  be  sinusoidally  distributed  as 

q  =  qQ  sin  ax  sin  By  (12) 

where  a  ^  n/a  ,  B  =  u/b. 

The  governing  equations  (7)  with  pressure  distribution  given  by  Eq. 

(12),  boundary  conditions  (9),  and  neutral -surface  location  criteria  (11)  are 
all  satisfied  exactly  in  closed  form  by  the  following  set  of  functions. 

u°  =  U  cos  ax  sin  By 

v°  =  V  sin  ax  cos  By 

w°  =  W  sin  ax  sin  ey  (13) 

i i  =  Y  sin  ax  cos  By 

iji  =  X  cos  ax  sin  By 

Then  differential  equation  set  (7)  reduces  to  algebraic  form  as 
follows: 

[Crs]{U,V,W,Y,X}T  =  (0,0,qo,0,0}T  (14) 

Here  [C  ]  is  a  symmetric  matrix  with  coefficients 

^11  =  A^a2  +  AggB2  ;  C12  =  (Ai2+Agg)aB  >  ^13  =  0  ; 

C 1 4  =  (B12+B66)aB  ;  C15  =  Bua2  +  B6632  ; 

C22  =  Agga2  +  A22S2  »  C23  =  0  ;  C24  =  Bgga2  +  B2262 

O5) 

C25  =  C14  ;  C33  =  K|s55a2  +  K4S44B2  ;  C34  =  K4S44B  ; 

C35  =  ^58550  ;  C44  =  Dgga2  +  D2262  +  K4S44  ; 

C45  *  (Di2+D66)a&  ;  C55  =  Dna2  +  D66a2  +  k|S55 


« 


Finite-Element  Analysis 
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Here  we  present  a  mixed  finite-element  model  associated  with  Eqs.  (1)- 
(  4  )  governing  the  bending  of  laminated,  thick  composite  plates.  The  word 
"mixed"  implies  that  independent  approximations  are  used  for  all  of  the 


variables,  u.v.w.ii'  >  and  ip  .  Using  the  thin-plate  equations  of  layered 
x  y 


composite  plates,  and  treating 


w,x  +  *x  =  0  ’  w’y  +  =  0 


(16) 


as  constraints,  Reddy17  presented  a  variational  formulation  of  Eqs.  (1)- 
(4).  That  is,  the  thick-plate  theory  can  be  interpreted  as  one  resulting 
from  the  thin-plate  theory  by  treating  the  slope-displacement  relations  as 
constraints.  The  Lagrange  multipliers  associated  with  these  constraints  are 


found  to  be  the  thickness-shear  stress  resultants,  Q  and  Q  .  The  model 

x  y 


described  here  is  essentially  the  same  as  in  Ref.  17. 

Suppose  that  the  region  occupied  by  the  plate  is  given  by  GX(-h/2,h/2) , 
where  ft  denotes  the  middle  plane  (x-y).  As  noted  earlier,  the  thickness 
direction  is  integrated  into  the  coefficients,  A..,  B..,  and  D.  ..  Hence,  we 

1  J  1  J  I  J 


divide  the  plate  into  a  finite  number  of  elements,  denoted  by  n  Over 


each  element  we  assume  that  the  variables  u,v,w,ijix,  and  ip  are  interpolated 


by  expressions  of  the  form 


n 


4>6  =  Z  N.  ^ 


(17) 


where  <j>e  denotes  the  restriction  of  a  typical  variable  to  ng,  its  value 


at  node  i  (of  element  n  ),  and  N .  are  the  linearly  independent  interpolation 


functions  associated  with  the  typical  element.  Since  we  are  concerned  here 
with  rectangular  plates,  the  typical  element  is  chosen  to  be  the  four-node 
(n=4)  quadrilateral  element  of  the  serendipity  family. 


.  MiBbltfltiii  I  , 


1 


—  - - 
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Substituting  expressions  of  the  form  (17)  for  u,v,w,^  ,  and  y  into 

*  y 

the  total  potential  energy  associated  with  the  case  of  a  cross-ply  laminate 

if  /_3y.\2  +  n  /iu>2  +  on  Miv  +  oft  iy.lv  +  A  (avx2 
2  J?  ^Au  ax'  A66(3y)  2A12  3x  3y  2A66  3y  3x  22l3y) 


.  a  /ll^2  +  l£/R  !!x  +  d  3^y\  +  d  /iy.  +  ivw^fx  +  ^y.\ 

A66(3y)  +  3x^Bu  3x  Bl2  3y  '  3y  3x^  3y  3x  ' 

+  IV/O  +  B  ffV  +  \B  au  +  B  av.\  +  d  #3i!x  +  +  i 1) 

+  3y^°12  3x  22  3y  3x  1 1 3x  123y^  Bs6*  3y  3x  3y  3x 

rs.i. 


+  ^lf^220)*Dn(^)a.  t,66(^,2 

a0x  3vv  dip  dip  2  3*v  2  /3w  ; 

+  2Di o — -  — ^  +  2Dfic  — -  — ^  +  022(— X-)  +  Dec ( — ^-)  +  A.c  (^  +  Ip  ) 
12  3x  3y  66  3y  3x  22V  3y  '  66V  3x  '  5:>v3x  x 

+  A44(|^  +  ■+  )2}  dxdy  +  qw  dx  dy  (18; 

we  obtain,  for  each  element, 

[Ke](Ae}  =  {Fe}  (19; 

where  Ue}  =  {u^v^.w®,^-  ,4'yi)T,  and 

Mj  ■  A>l^j  *  ^  -  ».2^  *  A66^ 

K!j  •  0,  =  Bn^j  ♦  Kj=  •  B12s‘*  ♦  B66S$ 

Kij  ■  *  A«*V  Kfj  ■  °-  K1j  ■  + 

Kt!  ■  B*4i +  B24p  Ku  5  Ass5u  +  A-*5*j>  Ku  *  *55S“ 

K?!  ■  A--Sir  Kil •  °>>su  *  *  A”s?j 

*  DseS$,  •  D^J  *  D24j  *  A„s“j 

F|  =  [  qN.  dx  dy,  F“  =  0,  a  =  1,2, 4, 5,  ^  =  KKg  Sofl,  (a, 6=4, 5) 


F?  =  q N.  dx  dy,  F“  =  0,  a  =  1,2, 4, 5,  A 

i  j  i 

e 

sfj  *  L  Ve  "j.n  d*  d*  (E.n-o,x,y),  ^  -  S 
Me 
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The  element  equations  (19)  are  assembled  in  the  usual  manner,  and 
boundary  conditions  are  applied  before  solving  the  equations. 

Numerical  Results 

As  the  first  example,  we  take  the  case  of  a  homogeneous  (single-layer) 
plate  of  transversely  isotropic  bimodulus  material.  The  plane  of  isotropy 
is  assumed  to  coincide  with  the  midplane  of  the  plate,  and  the  inplane 
Poisson's  ratio  is  assumed  to  be  zero.  Then  the  closed-form  solution  reduces 
to  the  simplified  form  presented  in  Appendix  B.  Numerical  results  are  pre¬ 
sented  in  Tables  1  and  2. 


Table  1.  Comparison  of  Neutral -Surface  Locations  for 
Transversely  Isotropic  Square  Plate 


Et'WGzc 

Neutral- 

-Surface  Location  Z 

W0-1 

0.3 

0.5 

Exact  Closed-Form  Solution: 

0.5 

-  0.08578 

-  0.08578 

-  0.08578 

1.0 

0 

0 

0 

2.0 

+  0.08578 

+  0.08578 

+  0.08578 

Simplified  Approximate  Solution: 

0.5 

-  0.08579 

-  0.08579 

-  0.08579 

1.0 

0 

0 

0 

2.0 

+  0.08579 

+  0.08579 

+  0.08579 

Mixed  Finite-Element  Solution: 

0.5 

-  0.08578 

-  0.08578 

-  0.08578 

1.0 

0 

0 

0 

2.0 

+  0.08578 

+  0.08578 

+  0.08578 
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Table  2.  Comparison  of  Maximum  Deflections  for  Transversely 
Isotropic  Square  Plate  (h/b=0.1,  K2i5/6) 


VEc=Gzt/6zc 

Dimensionless  Deflection  WEch3/qQb4 

W0-1 

0.3 

0.5 

Exact  Closed-Form  Solution: 

0.5 

0.05348 

0.04774 

0.04660 

1.0 

0.03688 

0.03283 

0.03201 

2.0 

0.02674 

0.02387 

0.02330 

Simplified  Approximate  Solution: 

0.5 

0.05004 

0.04660 

0.04591 

1.0 

0.03445 

0.03202 

0.03153 

2.0 

0.02530 

0.02342 

0.02296 

Mixed  Finite-Element  Solution: 

0.5 

0.05329 

0.04743 

0.04626 

1.0 

0.03675 

0.03261 

0.03178 

2.0 

0.02664 

0.02371 

0.02313 

It  is  noted  that  the  middle-surface  location  is  independent  of  Gzc  and 
Gzt>  The  agreement  among  the  results  obtained  by  all  three  solutions  is  quite 
good. 


As  examples  of  some  actual  bimodulus  materials,  aramid-cord/rubber  and 
polyester-cord/rubber  are  selected.  The  material  properties  used  are  listed  in 
Table  3.  The  data  are  based  on  test  results  of  Patel  et  al.3,  using  the  data- 


reduction  procedure  of  Model  2  in  Ref.  4,  except  for  the  thickness-shear  moduli, 


which  were  estimated  as  explained  in  Appendix  C. 


Table  3.  Elastic  Properties  for  Two  Tire-Cord/Rubber,  a 
Unidirectional,  Bimodulus  Composite  Materials 


Property  and  Units 

Aramid/Rubber 

k=l  k=2 

Polyester  Rubber 

k=l  k=2 

Longitudinal  Young's  modulus,  GPa 

3.58 

0.0120 

0.617 

0.0369 

Transverse  Young's  modulus,  GPa 

0.00909 

0.0120 

0.00800 

0.0106 

Major  Poisson's  ratio,  dimensionless^ 

0.416 

0.205 

0.475 

0.185 

Longitudinal-transverse  shear  modulus,  GPac 

0.00370 

0.00370 

0.00262 

0.00267 

Transverse- thickness  shear  modulus,  GPa 

0.00290 

0.00499 

0.00233 

0.00475 

aFiber-direction  tension  is  denoted  by  k=l ,  and  fiber-direction  compression  by 
h  k=2. 

It  is  assumed  that  the  minor  Poisson's  ratio  is  given  by  the  reciprocal 
relation. 

^It  is  assumed  that  the  longitudinal-thickness  shear  modulus  is  equal  to  this 
one. 

Numerical  results  for  single-layer  rectangular  plates  with  the  fibers 
oriented  parallel  to  the  x  axis  are  given  in  Table  4,  while  those  for  cross- 
ply  plates  (stacking  sequence  as  described  in  Appendix  A)  are  listed  in  Table 
5. 

As  can  be  seen  in  Tables  4  and  5,  the  agreement  between  the  closed-form 
and  finite-element  results  for  both  neutral -surface  position  and  deflection 
is  extremely  good.  Thus,  it  can  be  considered  that  the  finite-element  analysis 
has  been  soundly  validated,  and  can  now  be  used  for  more  complicated  combina¬ 
tions  of  loading,  geometry,  and  boundary  conditions  not  amenable  to  closed- 
form  solutions. 

It  is  noted  that  the  aramid-rubber  plates,  in  both  the  single-ply  and 
cross-ply  cases,  have  noticeably  larger  values  of  than  the  polyester- 
rubber  plates.  This  result  is  undoubtedly  due  to  the  more  pronounced  bi- 
modulus  effect  in  the  fiber-direction  Young's  modulus  of  the  aramid-rubber. 
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Table  4.  Neutral -Surface  Positions  and  Dimensionless  Deflections 
for  Rectangular  Plates  of  Single-Layer  0°  Aramid-Rubber 
and  Polyester-Rubber,  as  Determined  by  Two  Different 
Methods  (Thickness/Width,  h/b=0.1t;K2=5/6) 


Aspect 

Ratio 

Z 

X 

Z 

y 

1  ■■  ....  i.  .  .  .  _ 

WE22ch3/qob4 

C.F* 

F.E? 

C.F* 

F.E? 

C.F? 

F.E? 

Aramid-Rubber: 

0.5 

0.4453 

0.4454 

-  0.3304 

-  0.3007 

0.002544 

0.002750 

0.6 

0.4452 

0.4452 

-  0.2941 

-  0.2734 

0.004560 

0.004827 

0.7 

0.4447 

0.4447 

-  0.2564 

-  0.2419 

0.007393 

0.007712 

0.8 

0.4440 

0.4440 

-  0.2220 

-  0.2117 

0.01105 

0.01140 

0.9 

0.4431 

0.4431 

-  0.1923 

-  0.1846 

0.01545 

0.01582 

1.0 

0.4420 

0.4420 

-  0.1671 

-  0.1614 

0.02046 

0.02083 

1.2 

0.4394 

0.4394 

-  0.1285 

-  0.1250 

0.03160 

0.03193 

1.4 

0.4363 

0.4363 

-  0.1015 

-  0.09919 

0.04313 

0.04335 

1.6 

0.4328 

0.4329 

-  0.08228 

-  0.08070 

0.05406 

0.05416 

1.8 

0.4292 

0.4294 

-  0.06838 

-  0.06724 

0.06390 

0.06388 

2.0 

0.4253 

0.4254 

-  0.05813 

-  0.05727 

0.07250 

0.07236 

Polyester-Rubber: 

0.5 

0.3044 

0.3045 

-  0.1597 

-  0.1234 

0.001529 

0.001971 

0.6 

0.3044 

0.3045 

-  0.1538 

-  0.1245 

0.002652 

0.003265 

0.7 

0.3042 

0.3044 

-  0.1426 

-  0.1198 

0.004283 

0.005075 

0.8 

0.3039 

0  3041 

-  0.1299 

-  0.1124 

0.006517 

0.007487 

0.9 

0.3035 

0.3037 

-  0.1174 

-  0.1041 

0.009421 

0.01055 

1.0 

0.3029 

0.3031 

-  0.1061 

-  0.09586 

0.01303 

0.01430 

1.2 

0.3015 

0.3018 

-  0.08728 

-  0.08111 

0.02223 

0.02367 

1.4 

0.2999 

0.3001 

-  0.07329 

-  0.06941 

0.03348 

0.03492 

1.6 

0.2979 

0.2982 

-  0.06296 

-  0.06042 

0.04574 

0.04703 

1.8 

0.2957 

0.2960 

-  0.05528 

-  0.05356 

0.05793 

0.05897 

2.0 

0.2934 

0.2936 

-  0.04959 

-  0.04828 

0.06925 

0.07003 

★ 

C.F.  denotes  closed-form  solution;  F.E.  signifies  finite-element  solution. 
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Table  5. 

Neutral -Surface  Positions  and  Dimensionless  Deflections 
for  Rectangular  Plates  of  Two-Layer  Cross-Ply  Aramid- 
Rubber  and  Polyester-Rubber,  as  Determined  by  Two  Different 

Methods  (Thickness/Width,  h/b=0.1;  K2=5/6) 

Z 

X 

Z 

V 

WE22ch3/qob4 

Aspect 

Ratio 

C.  Ft 

F.E? 

C.F?  F.E? 

C.F? 

F.E? 

Aramid-Rubber: 

0.5 

0.4433 

0.4431 

-  0.06343  -  0.06223 

0.002472 

0.002576 

0.6 

0.4427 

0.4426 

-  0.05478  -  0.05443 

0.004388 

0.004518 

0.7 

0.4418 

0.4418 

-  0.04794  .  0.04778 

0.007072 

0.007220 

0.8 

0.4407 

0.4407 

-  0.04247  -  0.04237 

0.01054 

0.01070 

0.9 

0.4396 

0.4396 

-  0.03803  -  0.03795 

0.01475 

0.01490 

1.0 

0.4384 

0.4384 

-  0.03437  -  0.03430 

0.01957 

0.01972 

1.2 

0.4356 

0.4356 

-  0.02833  -  0.02860 

0.03043 

0.03054 

1.4 

0.4326 

0.4325 

-  0.02470  -  0.02477 

0.04185 

0.04190 

1.6 

0.4292 

0.4292 

-  0.02160  -  0.02165 

0.05282 

0.05280  j 

0.06264 

1.8 

0.4257 

0.4256 

-  0.01922  -  0.01923 

0.06277 

2.0 

0.4219 

0.4219 

-  0.01735  -  0.01734 

0.07151 

0.07137 

Polyester-Rubber: 

•j 

0.5 

0.3650 

0.3652 

-  0.1285  -  0.1256 

0.002539 

0.002732 

0.6 

0.3644 

0.3646 

-  0.1178  -  0.1164 

0.004527 

0.004772 

0.7 

0.3638 

0.3639 

-  0.1097  -  0.1089 

0.007288 

0.007575 

0.8 

0.3631 

0.3631 

-  0. 1036  -  0.1031 

0.01078 

0.01109 

0.9 

0.3622 

0.3622 

-  0.09886  -  0.09859 

0.01487 

0.01519 

1.0 

0.3613 

0.3613 

-  0.09526  -  0.09502 

0.01933 

0.01966 

1.2 

0.3593 

0.3593 

-  0.09001  -  0.09000 

0.02846 

0.02879 

1.4 

0.3571 

0.3570 

-  0.08660  -  0.08660 

0.03674 

0.03707 

1.6 

0.3546 

0.3545 

-  0.08430  -  0.08430 

0.04356 

0.04339 

1.8 

0.3519 

0.3518 

-  0.08267  -  0.08267 

0.04890 

0.04925 

2.0 

0.3491 

0.3490 

-  0.08150  -  0.08150 

0.05301 

0.05337 

*C.F.  denotes  closed-form  solution;  F.E.  signifies  finite-element 

solution. 
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Also,  it  is  interesting  to  observe  that  there  are  only  very  slight  differences 
in  Zx  and  deflection  in  going  from  a  single  layer  to  a  cross-ply  laminate. 

This  is  in  contrast  to  the  behavior  of  the  polyester  -  rubber  results  and  in 
considerable  contrast  to  ordinary  materials  (which,  of  course  have  Zx  =  0 
for  the  single-layer  case).  The  most  pronounced  change  in  going  from  the 
single-layer  to  the  two-layer  case  is  the  drastic  decrease  in  Zy  for  the 
a rami d- rubber. 

It  should  be  cautioned  that  in  the  case  of  the  closed-form  solution,  de¬ 
flections  due  to  various  sinusoidally  distributed  loadings  cannot  be  super¬ 
imposed  for  bimodulus-material  plates.  The  reason  superposition  is  not  valid 
here  is  that  the  necessary  conditions  for  homogeneity  of  neutral -surface  loca¬ 
tions  are  not  valid  under  superposition  conditions,  since,  in  general 


Z<’>  *!'>  (*,yMZ<2>  *<2>  (X,,) 


■7  _  n 

zn - TTT 


nr 

(x,y)  +  i|»XjX  (x,y) 


X  ,x 


f  constant 


even  though  zj^,  zj^,  ...  for  the  various  individual  Fourier  components  are 
constants.  However,  the  finite-element  solution  is  not  subject  to  these  limi¬ 
tations,  since  it  provides  for  stepwise  variation  in  neutral -surface  location. 
Fig.  1  shows  results  for  both  sinusoidally  and  uniformly  distributed  loadings. 

Concluding  Remarks 

Both  finite-element  and  closed-form  solutions  have  been  found  for 
thick,  rectangular  plates  of  single-layer  and  cross-ply  laminates  of  bi¬ 
modulus  materials.  Excellent  agreement  was  obtained,  and  thus  the  finite- 
element  formulation  of  this  problem  is  considered  to  have  been  validated 
against  an  accurate  benchmark. 
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Effect  of  plate  aspect  ratio  on  dimensionless  displacements  of  aranid 

rubber  and  polyester-rubber  rectangular  plates  under  uniform  and  sinu 
soidally  distributed  loadings  (UL  and  SSL).  Here  w=wE-i  h5/q  b-  and 
u=103U/W.  “c  0 


The  research  reported  here  is  currently  being  extended  to  (1)  thermal 
bending  due  to  changes  in  midplane  temperature  and  in  gradient  through  the 
thickness,  and  (2)  free  vibration. 
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Appendix  A:  Derivation  of  the  Plate  Stiffnesses  for  a  Two- 
Layer  Cross-Ply  Laminate  of  Bimodulus  Material 

In  laminates  containing  bimodulus  materials,  the  results  of  evaluating 
the  integrals  for  the  plate  stiffnesses,  Eqs.  (5),  are  more  complicated  than 
those  for  ordinary-material  laminates,  since  the  individual-layer  plane- 
stress-reduced  stiffnesses  depend  upon  the  neutral -surface  location.  Here 
we  derive  the  expressions  for  the  case  of  a  two-layer,  cross-ply  laminate. 

Each  layer  has  the  same  thickness  and  the  same  bimodulus  orthotropic  elastic 
properties  with  respect  to  the  fiber  direction.  The  bottom  layer  is  denoted 
as  layer  1,  i.e.,  I  =  1  in  is  oriented  in  the  x  direction,  and  occupies 

the  thickness-direction  interval  from  z  =  0  to  z  =  h/2,  where  z  is  measured 
position  downward  from  the  midplane.  The  top  layer  (i= 2)  is  oriented  in  the 
y  direction  and  is  located  from  z  =  -  h/2  to  z  =  0.  In  the  case  considered, 
it  is  assumed  that  the  upper  portion  of  the  top  layer  (i=2)  is  in  compression 
(k=2  in  Qijk£)  in  the  fiber  direction  and  that  the  lower  portion  of  the  top 
layer  is  in  tension  (k=l),  while  the  inner  portion  of  the  bottom  layer  U=l) 
is  in  compression  (k=2)  and  the  outer  portion  of  this  layer  in  tension  (k=l). 
Summarizing,  the  four  regions  are  as  follows: 


Layer 

Reqion 

Fiber-Direction 
Tension  or  Compression 

i  =  2 

-  h/2  to  t„y 

Compression  (k=2) 

i  -  2 

t0  0 

Tension  (k=l) 

i  =  1 

0  t0  2nx 

Compression  (k=2) 

t  =  1 

2nx  t0  W 

Tension  (k=l) 

It  is  noted  that  it  is  assumed  that  the  x-direction  neutral -surface  location 
znx  i  0,  while  the  y-direction  one  (zny)  is  negative. 


i 


21 


Thus,  the  integral  for  A.,  in  Eqs.  (5)  is  subdivided  into  four  regions, 

'  J 

in  each  of  which  the  plane-stress  reduced  stiffnesses  are  constant. 

,2ny  r°  rZnx  ,h/2 

A. .  =  Q. .  dz  +  Q.  .  dz  +  Q. .  dz  +  0. .  dz 

1J  J-h/2  1J22  Jz  1J12  JO  1J21  K  ljn 

ny  nx 

=  (Q. .  +  Q . .  ) (h/2)  +  (Q . .  -Q..  )z  +  (Q.  .  -Q..  )z  (4 

vwijn  yij22,v  '  vyij2i  yijny  nx  v w i j 2 2  wiji2;  ny 


Introducing  Z  =  z  / h  and  Z  =  z  / h,  one  obtains 
x  nx  y  ny 


Aij/h  -  (l/2)(Q.j j22  +  +  (Qij21-  QijnJZx  +  ^ij22  ‘  ^iji2^Z- 


In  similar  fashion,  the  next  two  integrals  in  Eqs.  (5)  become 


4Bij/h2 


(1/2)(Q- 


J22  ^ 


2  ( Q .  . 

VX1J21 


-  Q-  •  )Z2 
yijn'  x 


2(Q. . 
1J22 


-  Q. .  )Z2 

viji2y  y 


12D.  ./h3  *  ( 1/2 XQ -  •  +Q--  )  +  4(0.  .  -Q..  )Z3  +  4(0.  .  -Q..  )Z3 

1J  '  A  IJll  y1J22'  V J2 1  X  VW1J21  yl J 12 '  y 


The  expression  for  S .  ./h  is  the  same  as  for  A../h,  Eq.  (A2). 

'  J  1  J 

To  apply  Eqs.  (A2)-(A4)  to  a  single-layer  plate  with  the  fibers  oriented 

in  the  x  direction,  it  is  necessary  to  merely  set  Z^  =  0.  In  deriving  these 

equations,  it  was  assumed  that  Z  >  0  and  Z  <_  0.  In  the  event  that  the  final 

«  y 

results  obtained  for  Z  and  Z  did  not  meet  these  conditions,  obviously  Eqs. 

x  y 

(A2)-(A4)  would  not  be  valid.and  it  would  become  necessary  to  investigate 
another  of  the  other  three  possible  cases: 


Z  >  0  and  Z  >  0 
x  y 

Z  <  0  and  Z  <  0 
x  y 

Z  <  0  and  Z  >  0 
x  y 
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Fortunately,  however,  in  all  of  the  static  problems  treated  here,  the 
conditions  for  the  case  derived  in  this  appendix  are  satisfied. 

Appendix  B:  Reduction  of  Closed-Form  Solution  for  Single- 
Layer  Transversely  Isotropic  Material 


In  order  to  obtain  simple,  concise  expressions  for  the  neutral-surface 
location  and  the  deflection,  the  closed- form  equations  given  in  the  body  of 
the  paper  are  reduced  to  the  special  case  of  a  square  plate  made  of  a  trans¬ 
versely  isotropic  bimodulus  material  with  an  inplane  Poisson's  ratio  of  zero. 

Then 


An  =  A22 

=  A, 

Ai2  =  0, 

^66  =  A/2 

Bn  =  B22 

=  B, 

B 1 2  =  0, 

B66  =  B/2 

Dll  =  °22 

=  D, 

0i2  =  0, 

066  =  D/2 

S44  =  S55 

=  s. 

S4  5  =  0 

(81) 


Now  Eqs.  (15)  reduce  to  the  following,  since  B  =  a: 


Cn  =  (3/2)Aa2, 

C12  *  (l/2)Aa2 , 

Cis  =  0 

C14  =  ( l/2)Ba2 , 

Ci 5  =  (3/2)Ba2 , 

Cji  =  Cij’ 

C22  =  Cn»  c23 

=  0,  C24  =  Ci 5 , 

C25  =  Ci4, 

(B2) 

C33  =  2K2Sa2 ,  C34  =  K2$a,  C35  =  C34, 

C44  =  (3/2)Da2  +  K2S,  C45  =  (l/2)Da2,  C55  =  C44 

The  biaxial  symmetry  of  this  special  case  requires 

b  s  a,  8  -  a,  V  -  U.  Y  -  X.  2nx  =  2ny  =  2n 


(B3) 


Using  Eqs.  (B2)  and  (B3)  in  the  first  two  of  Eqs.  (14),  one  finds  that 
for  this  special  case 


nx 


■ B/A 


(B4) 


4  "  i'”" 
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Using  the  fourth  equation  of  Eqs.  (14),  one  obtains 

X/"  ■  *  K*S  <B5> 

It  is  noted  that  the  bending  slope  vanishes  for  both  S  =  0  and  1/S  =  0. 
Finally,  the  third  of  Eqs.  (14)  yields 

W  =  4(D-BZn)a4  ^  +  2( D-Bz^) (cx2/K2s ) ]  (B6) 

It  is  seen  that  the  quantity  in  front  of  the  first  bracket  on  the  right  side 
of  Eq.  (B6)  is  equal  to  the  deflection  of  a  thin  isotropic  plate.  The  second 
term  inside  the  brackets  is  the  fractional  increase  in  deflection  due  to 
thickness-shear  deformation,  which  obviously  increases  as  Gz  is  decreased. 

The  quantity  D-Bzn  is  the  so-called  reduced  stiffness,  first  obtained  for 
laminated,  isotropic  thin  plates  by  Pister15. 

For  the  present  case,  Eqs.  (A2)-(A4)  become 


A/h  =  Q  +  ZaQ  ,  4B/h2  =  -  (1/2)aQ  (1-4Z2) 

12D/h3  =  Q  +  4Z3AQ 


(B7) 


Here 

Q  =  (1/2)(QC+Qt)  ,  AQ  =  Qc  -  Qt  (B8) 

Combining  Eqs.  (B4)  and  (B7),  one  obtains  the  following  quadratic 
expression  for  Z: 

Z  =  -(Q/AQ)  +  [(Q/aQ)2  -  (1/4)]*  (B9) 

Also 

S/h  =  (V2)(Gzc+Gzt)  +  Z(Gzc-Gzt)  (BIO) 


i 
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Appendix  C:  Method  of  Estimating  Transverse-Thickness 
Shear  Moduli 


In  the  tests  reported  by  Patel  et  al.3, only  inplane  compliances  were 
measured.  Thus,  it  is  necessary  to  estimate  the  values  for  the  thickness- 
shear  moduli,  which  are  needed  for  the  thick-plate  analysis. 

It  is  believed  to  be  a  reasonable  engineering  assumption  to  assume 
that  an  individual  composite-material  layer  is  transversely  isotropic  with 
the  plane  of  isotropy  being  the  cross-sectional  plane,  i.e.,  the  plane  normal 
to  the  fibers.  Thus,  it  follows  that  the  longitudinal-thickness  shear  modulus 
(G13)  is  equal  to  the  inplane  (longitudinal-transverse)  shear  modulus  (Gn). 

Estimation  of  the  other  thickness-shear  modulus  is  more  complicated.  One 
can  use  the  well-known  isotropic  relation  for  the  transverse-thickness  shear 
modulus  G2 3  in  terms  of  the  thickness  Voung's  modulus  E3  and  transverse¬ 
thickness  Poisson's  ratio  v23  provided  the  latter  two  quantities  are  known: 

G23k  ■  V^^k”  ‘C’> 


In  view  of  the  transverse-isotropy  assumption  mentioned  above,  it  can 
be  assumed  that  E3k  =  E2k,  which  was  obtained  from  the  inplane  tests. 

Foye16  presented  a  relation  for  v23  which  can  be  rewritten  in  the 
following  form,  which  is  more  convenient  for  the  present  purpose: 


v  .  + 
12k 


[(vm/Em)  ~  (vf/Ef)]  vm(1~vm) 

l-vf,  ,  V_  Vf 


Ef  vf 


m  f 
Em  Ef 


It  is  noted  that  Ref.  3  provided  data  for  Ef,  Em,  and  V^.  It  is 

reasonable  to  use  a  value  of  0.499  for  v  of  rubber.  Thus,  the  only  unknown 

m 

quantity  on  the  right  side  of  Eq.  (C2)  is  which  could  be  computed  from 


25 


the  following  rule-of-mixtures  expression  for  v12^,  which  is  known  to  be 
very  accurate  for  polymer-matrix  composites: 

vi2k  =  vfkVf  +  vm(1-Vf}  (C3: 

The  rule-of-mixtures  expression  for  the  longitudinal  Young's  modulus 
is  also  known  to  be  accurate  for  polymer-matrix  composites: 

EU  ■  Efkvf  *  v-v  <c«: 

Unfortunately ,  however,  for  the  data  of  Ref.  3,  the  measured  values 
of  Elt  (tension)  were  higher  than  predicted  by  Eq.  (C4).  Thus,  it  was  de¬ 
cided  to  use  Eq.  (C4)  to  obtain  an  effective  fiber  volume  fraction  V^.,  and 

then  to  use  this  effective  value  to  predict  Efc  ( compress ion)  and  and 

★ 

Vfc-  However,  using  either  or  in  Eq.  (C3),  one  obtains  negative  vft 
and  values,  which  are  not  reasonable  physically.  Thus,  it  was  assumed 
that  due  to  the  loose  nature  of  the  cord,  that  it  was  not  restrained  by  the 
matrix.  Thus,  instead  of  obtaining  from  Eq.  (C3),  it  was  obtained  from 

vfk  ■  Vi2k/'vf  <C5) 

Sample  calculations  for  aramid-rubber  in  compression  are  as  follows. 
From  Eq.  (C4)  for  k  =  t: 


Vf  =  (Eit'Em^/(Eft"Em)  =  (3-58  *  °-008°)/(24-8  '  0.0030)  =  0.144 


Then,  using  Eq.  (C4)  for  k  =  c: 

Efc  *  -Eic"  Em(1-Vf^/Vf  =  C°-00120  '  0. 0080(0. 856)]/0. 144 
=  0.0358  GPa 
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From  Eq.  ( C 5 )  with  k  =  c, 

vfc  ■  Vc/vJ  =  0-205/0.144  =  1.42 
Using  Eq.  (C2),  one  obtains  v,3c  *  0.202 
Finally,  from  Eq.  (Cl).  G ,  =  0.00499  GPa 

.  )  L 
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